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^ . ABSTRACT 

^ ' We present chemical abundance measurements from medium resolution observations of 8 sub- 

damped Lyman-a absorber and 2 strong Lyman-limit systems at z < 1.5 observed with the 
MIKE spectrograph on the 6.5m Magellan II Clay telescope. These observations were taken 
as part of an ongoing project to determine abundances in Zabs"^ 1-5 quasar absorption line 
systems (QSOALS) focusing on sub-DLA systems. These observations increase the sample 
of Zn measurements in Zabs"^ 1-5 sub-DLAs by ^50%. Lines of Mg I, Mg II, Al II, Al III, 
Ca II, Mn II, Fe II, and Zn II were detected and column densities were determined. Zn II, 
a relatively undepleted element and tracer of the gas phase metallicity is detected in two of 
these systems, with [Zn/H]= — 0.05±0.12 and [Zn/H]>H-0.86. The latter system is however 
a weak system with Nh i <18.8, and therefore may need significant ionisation corrections to 
the abundances. Fe II lines were detected in all systems, with an average Fe abundance of 
<[Fe/H]>= —0.68, higher than typical Fe abundances for DLA systems at these redshifts. 
This high mean [Fe/H] could be due to less depletion of Fe onto dust grains, or to higher 
abundances in these systems. We also discuss the relative abundances in these absorbers. The 
k> ■ systems with high metallicity show high ratios of [Mn/Fe] and [Zn/Fe], as seen previously in 

r^ , another sub-DLA. These higher values of [Mn/Fe] could be a result of heavy depletion of Fe 

■ onto grains, unmixed gas, or an intrinsically non-solar abundance pattern. Based on Cloudy 

modeling, we do not expect ionisation effects to cause this phenomenon. 
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1 INTRODUCTION 

There are still many open questions concerning the processes of 
galaxy formation and evolution. Measurements of the abundances 
of heavy elements in galaxies give important information about 
the ongoing processes of star formation and death and the overall 
chemical enrichment of these galaxies. Studying galaxies at higher 
redshifts through emission is often difficult and time consum- 
ing, and biases towards more luminous galaxies are often present. 
Quasar absorption line systems (QSOALS) provide a means of 
studying the interstellar medium (ISM) of high redshift galaxies in- 
dependent of the morphology of galaxies (that is, not based on the 
selection of galaxies of certain morphology). In addition, absorp- 
tion lines in QSO spectra allow us t o study the diffuse in tergalactic 
medium using lines of O VI (e.g.. ISimcoe et alj|2002l) . or using 



X-Ra y absorption lines of more highly ionized species jpang et al] 
l2002h . 

Quasar absorption line systems with strong Lyman-a lines are 
often divided into two classes: Damped Lyman-Q (DLAs, log Nh i 
^ 20.3) and sub-Damped Lyman-a (sub-DLA 19 < log Nh i < 
20.3, Peroux et al. 2001) are expected to contain a major fraction 
of the neutral gas in the Universe, while the majority of the baryons 
are thought to lie in the highly ionized and diffuse Lyman-Q forest 
clouds with log Nh i < 14 in intergalactic space jPetitiean et al.l 
[19930. The DLA and sub-DLA systems are generally believed to 
be associated directly with galaxies, at all redshifts in which they 
are seen. Observations of the galaxies in emission responsible for 
the absorption line systems have been met with mixed results. In 
several instances, a galaxy at the same redshift in emission could 
not be found near the line of sight of the QSO (see for exam- 
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plelKulkami et alj|200d. I2001I : Ichen and Lanzett j|2003l : iRao et alj 



l200i ) 

DLAs have been the preferred systems for chemical abun- 
dance studies thanks to their high gas content. However, most 
DLAs are found to be metal po or, typically far below the solar 
level (e.g., iKulkami et iZI ( l2005h and references there in). Based 
on Fe II abundance measurements, |Perouxetal]j2003ah suggested 
stronger metallicity evolution in s ub-DLA systems than DLA sys- 
tems. Specifically, IPeroux et al. I j2003ah showed that the Nh i- 
weighted mean Fe metallicity in sub-DLAs increased from ^ 1/100 
Zq at z~4.5 to ~ 1/3 Zq at z~0.5. This has been validated 
with Z n, the more reliable metallicity inidicator by Kul karni et al.l 
( |2007|) . 

To date, few observations have been made of z < 1.5 sub- 
DLAs due to the lack of spectrographs with enough sensitivity at 
short wavelengths and the the paucity of known sub-DLAs in this 
redshift range. Redshifts z < 1.5 span 70% of the age of the Uni- 
verse (assuming a concordance cosmology of Jim = 0-3, Qa = 
0.7). This redshift range is clearly important for understanding 
the nature of sub-DLA systems and galactic chemical evolution as 
well. 

Recently IPer oux et al.l (12003 j. l2006al) : lYork et al.l 
Kulkarni et al. (2Q0'^: |Prochaska et alJj2'ooa) : lKhare et alj 



200e); 



20071) : 



Meiring et al., 1 20071) have suggested that sub-DLAs may be more 



metal rich than DLAs, and may contribute significantly to the cos- 
mic metal budget. Several sub-D LA systems have actually been 



seen w i th super-solar abundan c es iKhare et al. 



l2006al : iProchaska et"al] I2OO6I : iMeiring et all 



2004l ; |Perouxetal.l 
20071) . Super-solar 



abundances have also been reported in strong Lyman-a forest lines 
with log Nh i <16, although ionisati on correction factors are im- 
portant and not easily constrain ed ( Charlt on et alj2003l : lDing et al.l 
l2003l : lMasiero et ali20o3 : [Schave 2 007). 

As the DLA systems dominate the Nh i weighted mean metal- 
licity, previous studies have been dominated by DLAs. At least 
an equally large sample of sub-DLAs is needed to determine the 
overall metallicity evolution of QSO absorbers. The H I column 
density distribution for QSO absorbers shows that systems with 19 
< log Nh I < 20.3 are ~8 times more numerous than the higher 
column density classical damped Lya systems, making sub-DLAs 
more readily available probes of neutral gas in the distant Universe. 

Many elements have been detected in QSO absorber systems, 
including C, N, O, Mg, Si, S, Ca, Ti, Cr, Mn, Fe, Ni, and Zn. Zn is 
often the preferred tracer of the gas-phase metallicity as it is rela- 
tively undepleted in the Galactic ISM, especially when the fraction 
of H in molecular form is low, as is the case in mo st DLAs. Zn 
also tr acks the Fe abundance in Galactic stars (e.g., iNissen et al.l 
( l2004 l)). and the lines of Zn II AA 2026,2062 are relatively weak 
and typically unsaturated. These lines can be covered with ground 
based spectroscopy over a wide range of redshifts, 0.65 < z < 3.5, 
which covers a large portion (~ 45%) of the history of the Uni- 
verse. Abundances of refractory elements such as Cr and Fe rel- 
ative to Zn also give a measure of the amount of dust extinction 
jYork et alj|2006l) . 

The number of sub-DLAs studie d to date is still s mall, 
with the largest sam ples coming from iPeroux et al.l ( l2003al) and 
iMeiring et al.l j2007l) . In this paper we add 10 new systems to 
our previous medium-resolution measurements of sub-DLAs taken 
with the Magellan Inamori Kyocera Echelle (MIKE) spectrograph 
on the 6.5m Clay telescope at Las Campanas observatory. As the 
amount of published data on abundances of low redshift sub-DLAs 
is still small, we present these observations here with a more de- 
tailed analysis of the abundances (both absolute and relative) and 



kinematics in a future paper with our complete sample. The struc- 
ture of this paper is as follows. In § 2, we discuss details of our 
observations and data reduction techniques. In § 3, column density 
determinations are discussed. § 4 gives information on the individ- 
ual objects, and in § 5 we discuss the abundances of these absorbers. 
We also provide an appendix at the end of this paper, showing plots 
of the UV spectra with the fits to the Lyman-a lines. 



2 OBSERVATIONS AND DATA REDUCTION 

The observations presented here were made with the 6.5m 
Magellan-II Clay telescope and the Magellan Inamor i Kyocera 
Echelle (MIKE) spectrograph teemstein et"d] l2003h in 2007 
March. This is a double sided spectrograph with both a blue and a 
red camera, providing for simultaneous wavelength coverage from 
~3340 A to ~9400 A. Targets were observed in multiple exposures 
of 1800 to 2700 sec each to minimize cosmic ray defects. The see- 
ing was typically < I", averaging ~ 0.7". All of the target QSOs 
were observed with the I"x5" slit and the spectra were binned 
2x3 (spatial by spectral) during readout. The resolving power of 
the MIKE spectrograph is ~ 19,000 and ~25,000 on the red and 
blue sides respectively with a 1" slit after binning. Table 1 gives a 
summary of the observations. 

These spectra were reduced using the MIKE pipeline reduc- 
tion code in IDL developed by S. Buries, J. X. Prochaska, and R. 
Bernstein. Wavelengths were calibrated using a Th-Ar comparison 
lamp taken after each exposure. The data were first bias subtracted 
from the overscan region and flat-fielded. The data were then sky- 
subtracted and the spectral orders were extracted using the traces 
from flat field images. These extracted spectra were then corrected 
for heliocentric velocities and converted to vacuum wavelengths. 
Each individual order was then extracted from the IDL structure 
created by the pipeline software and combined in IRAF using re- 
jection parameters to reduce the effects of cosmic rays. These com- 
bined spectra were then normalized using a polynomial function to 
fit the continuum. These functions were typically of order five or 
less. 

Our sample consists of 10 absorbers, 8 with 19.0 < log Nh i 
< 20.3 and 2 absorbers with log Nh i < 19.0, at Zabs^ 1-5, where 
few observations have been made to date. We selected these sys- 
tems based on several criteria. Firstly, the absorbers were chosen to 
have 0.65 < Zabs^ L5. Faint QSOs (mv < 18.8) were avoided. 
All absorbers have known Nh i, with 19.0<Nh i <20.3. Systems 
that were previously known to have extremely weak lines of Mg 
II (Wl''^^ < 0.6 A), which are very unlikely to hold detectable 
Zn II lines at the thresholds available to this program were also 
excluded. Systems with low Nh i (Nh i < 19.2) and faint QSOs 
(mv > 18.5) were also avoided. Throughout this paper the QSO 
names are given in J2000 coordinates, except in table 1 where the 
original name, based on J 1 950 coordinates is given as the name 
( iHewitt fc Burbidgelll987h . 

We have attempted to be as unprejudiced as possible when 
selecting targets for ob s ervatio n. Of the ~45 absorbers from 
iRao. Turnshek, fc Nestod ilOOdt) that are accessible with MIKE 
(i.e. had bright enough background QSOs and were observable 
within the RA and declination constraints), we have observed 38 
of these absorbers to date (Meiring et al. (2006, 2007), this work, 
and recent unpublished observations). Of the 7 that have not been 
observed, one was not selected because of high (yet technically 
possible) declination, one was excluded due to very low absorp- 
tion redshift which would make detection of the Zn II lines dif- 
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ficult at the wavelength cutoff of MIKE, and most of the others 
were skipped due to scheduhng and observing constraints. Only 
two of these 7 absorbers have W2796 < 1.0 A, and none have 
W27g6 < 0.5 A (i.e. the remaining systems are fairly strong Mg II 
syste ms), showing that this is a relatively unbiased sample. Further- 
more, lKulkamietai] ( l2007 l) showed that systems with large W2796 
have a wide range of [Zn/H], and that systems with high [Zn/H] 
have a large range of W2796, or that selection by W2796 in the low 
z Nh I surveys does not appear to significantly bias the sample to- 
wards higher metallicity systems. 

3 DETERMINATION OF COLUMN DENSITIES 

Column densities were determined from profile fitting with the 
package FITS6P tWelty et al. 199 1), which has evolved from the 
code by IVidal -Madiar et alj ( Il977h . FITS6P iteratively minimizes 
the value between the data and a theoretical Voigt profile that 
is convolved with the instrumental profile. The profile fit used mul- 
tiple components, tailored to the individual system. For the cen- 
tral, core components, the effective Doppler parameters (6e//) and 
radial velocities were determined from the weak and unsaturated 
lines, typically the Mg I A 2852 line. For the weaker components 
at higher radial velocities the b^ff and component velocity val- 
ues were determined from stronger transitions such as the Fe II 
AA 2344, 2382 lines and the Mg II AA 2796, 2803 lines. A set of 
6e// and v values were thus determined that reasonably fit all of 
the lines observed in the system. The same b^f / values were used 
for all the species. The atomic data used in line identification and 
profile fitting are from lMortonI ( |2003|) . 

In general, if a multiplet was observed, the lines were fit si- 
multaneously until convergence. For all of the systems, the Fe II A 
2344, 2374, 2382 lines were fit simultaneously to determine a set of 
column densities that fit the spectra reasonably well. Similarly, the 
Mg II AA 2796, 2803 lines were also fit together. Significant satu- 
ration of the Mg II A A 2796, 2803 and Al II A 1670 lines allowed 
for only lower limits to be placed on the column densities for these 
species. The Zn II A 2026.137 line is blended with a line of Mg I 
A 2026.477. The Mg I contribution to the line was estimated using 
the Mg I A 2852 line, for which fA is ~32 times that of the Mg I A 
2026 line. The Zn II components were then allowed to vary while 
the Mg I components were held fixed. Ncr 11 was determined by 
simultaneously fitting the Cr II A 2056 line and the blended Cr II -1- 
Zn II A 2062 line, where the contribution fr om Zn II was estim ated 
from the Zn II + Mg I A 2026 line. See also lKhare et"ai] ( |2004|) and 
iMeiring et al.l j2007l) for a discussion of the profile fitting scheme. 
In this paper we adopt the standard notation: 

[X/Y] = log(Nx/NH i) - log(X/H)0 (I) 

Solar system abundances have been adopted from lLodderj ( l2003h . 

We present the rest-frame equivalent widths (Wo) of the lines 
in Table 2. The la errors for the equivalent widths are given also 
and reflect both uncertainties in the continuum level and in the pho- 
ton noise. If a certain line was not detected, the limiting equivalent 
width was determined from the local signal to noise ratio (S/N) 
and uncertainties in the continuum placement based on a 3-pixel 
element. Assuming a linear curve of growth, a corresponding 3a 
upper limit was also placed on the column density. 

We give total column densities (the sum of the column den- 
sities in the individual components determined via profile fitting 
method) for these systems in table 3. Column densities were also 



obtained via the apparent optical depth (AOD) method as a consis- 
tency check, and are also listed in table 3. F or a discussion of the 
AOD method, see ISavagefcSembach|jl996 '). Errors for the AOD 
column densities include contributions from both the photon noise 
and continuum placement. In most cases, the column densities from 
the profile fitting and AOD methods agree to within the error bars, 
especially for the weak and unsaturated lines. 



4 DISCUSSION OF INDIVIDUAL OBJECTS 
4.1 Q1037+0028 (z„„ = 1.733) 

(System A:Zni,s=l-4 244): This QSO harbors a sub- DLA with log 
Nh I = 20.04±0.I2 jRao. Tumshek. fc Nestoill2006l) . We detected 
absorption lines for Mg II, Al II, Al III, Si II, Mn II and Fe II in this 
sub-DLA system. The Mg I A 2852 line was blended with telluric 
features at ~6916 A. A total of II components were used in the 
profile fitting for this system. No Zn II lines were detected in this 
system at S/N '--^ 16 in the region, and a limiting EW of Wo < 21 
mA. There were two small features that lined up with the compo- 
nents at 1 1 and 44 km s~^ from other lines, however, the equivalent 
width of the features was significant at only a <2a level. Total col- 
umn densities for this and all other systems are given in Table 3. 
The abundance based on Fe lines was [Fe/H]= — 0.67±0. 12. This 
system shows Al III/Al II < —0.96, a typical level of ionisation for 
a sub- DLA (e.g.. lDessauges-Zavadskv et alj l [2003h : lMeiring et al.l 
( I2OO7I) ). Rest frame equivalent widths for the lines in this system 
are found in table 2, along with the results of the profile fitting in ta- 
ble 4. Velocity plots of several lines of interest are shown in Figure 
I. 

There also appears to be a weaker system at Zats =0.6572, with 
Wo(Mg II 2796)=354 mA. We detect the Fe II AA 2344, 2382, 
2586, and 2600 lines along with the Mg I A 2852 line as well. The 
Fe II A 2374 line is blended with the Galactic Ca II A 3934 line 
for this system. Nh i cannot be determined for this weaker system 
due to the break in the continuum from the sub-DLA system at 
Zabs=1.4244 beyond the Lyman limit, so it is not included in our 
sample (which depends on known Nh i, as noted earlier). 



4.2 Q1054-0020 (z^^ =1.021) 

(System A:Za6s=0.8301): S ystem A is a strong Lyman-lim it sys- 
tem with Nh I =18.95±0.18 jRao. Tumshek. fc Nestoj2006h . Five 
components were used in the profile fits for system A. Lines of Mg 
I , Mg II, Mn II, and Fe II were all detected. No Zn II lines were 
seen at S/N^23 in the region. We did not detect Zn II A 2026 with 
a limiting EW of Wo< 10 mA, and [Zn/H]< +0.18. This system 
has a ratio of Nfc ii/Nh i= — 0.09±0. 18 dex. The true ratio may 
be lower due to ionisation (see our later comments on ionisation 
corrections at low Nh i). We also detected lines of Mn II AA 2576, 
2594 with [Mn/Fe]= -0.09±0.07. The Al II A 1670 line lies below 
the wavelengths covered by MIKE, so the Al III/Al II ratio could 
not be measured in this system. We did however detect strong Al 
III AA 1854, 162 lines. Velocity plots of several lines of interest are 
shown in Figure 2. Rest frame equivalent widths for the lines in this 
system can be found in table 2, with the results of the profile fitting 
analysis in table 5. 

(System B:Za(,s=0.9514): A total o f 7 components were used 
to fit system B with Nh i =19.28±0.02 jRao, Tumshek. fc Nesto3 
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Table 1. Summary of observed targets and their absorption systems. 



QSO 


Alternate name 
or SDSS ID 


mv or g 


RA 


Dec 




'^abs 


log Nh I 


E(B— V)g_i'' Exposure Reference 
time 


J2000 






J20()0 


J2000 






cm 


s 



Q1037+0028 


SDSS J103744.44+002809.2 


18.4 


10:37:44.40 


-1-00:28:09.2 


1.733 


1.4244 


20.04±0.12 


0.04 


10800 


1 


Q1054-0020 


SDSS J105440.98-002048.4 


18.3 


10:54:40.98 


-00:20:48.4 


1.021 


0.8301 


18.95±0.18 


. . i 


10800 


1 














0.9514 


19.28±0.02 


0.05 






Q1215-0034 


Q1213-002 


17.0 


12:15:49.81 


-00:34:32.1 


2.691 


1.5543 


19.56±0.02 


0.02 


5400 


2 


Q 1220-0040 


SDSS J122037.00-004032.4 


18.5 


12:20:37.00 


-00:40:32.4 


1.411 


0.9746 


20.20±0.07 


-0.02 


8100 


1 


Q1228-I-1018 


Q 1226+ 105 


18.5 


12:28:36.80 


-1-10:18:41.7 


2.305 


0.9376 


19.41±0.15 


0.00 


7200 


3 


Q1330-2056 


Q 1327-206 


17.0 


13:30:07.70 


-20:56:16.4 


1.165 


0.8526 


19.40±0.02 


N/A'= 


3600 


4,5 


Q1436-0051 


SDSS J143645.06-005150.6 


18.5 


14:36:45.03 


-00:51:50.6 


1.275 


0.7377 


20.08±0.11 


0.08 


5400 


1 














0.9281 


<18.8 


. . b 






Q1455-0045 


SDSS J145508.14-004507.4 


18.0 


14:55:08.14 


-00:45:07.5 


1.378 


1.0929 


20.08±0.06 


0.02 


8100 


1 



REFERENCES: (1) [Ne stor (20oj); (2) ILanzetta. Tumshek. & 
iBergeron. Kunth. fc mjdoricd < 198?!) 

"Determined using t he colour excess A{g — i) = {g — i) - 
[Richards etaP <200lh and the relation E{B -V)g-i = A{g - 
'The extintinction is attributed to the system with the higher Nh 
'^This QSO was not observed in the SDSS 



: Wolfd h987l) : (3) iBarthel. Tvtler. fc Thompsor] il990l) : (4) iKunth fc BergerorJ il984h : (5) 



- {g — i)med^ where {g — i)med is the median QSO colou r at the emission redshift from 
0(1 + ^a6s)"^-^/1.506 iKhare et alj2004l : lYork et al.l2006ll . 
J, although the true extincition could be due to either, or both systems. 



I2OO6I) . Only lines of Mg I, Mg II, and Fe II were detected in this sys- 
tem. We were able to constrain the metallicity of this system based 
on the non detection of the Zn II A 2026 line to be [Zn/H]< —0.21 
with S/N~30 in the region. This system also has a low metallicity 
based on Fe with [Fe/H]= — 1.09±0.02, although dust depletion 
may be important. The Al II A 1670 line lies below the covered 
wavelengths. The Al III A A 1854, 1862 lines were covered but not 
detected with Wo<22 mA. Velocity plots of several lines of interest 
are shown in Figure 3. The results of the profile fitting analysis are 
found in table 6, and the rest frame equivalent widths for the lines 
in this system in table 2. 



were used to fit the line profiles for this system. The Al II A 1670 
line lines below the wavelengths accessible to MIKE, so we were 
unable to determine the Al III/Al II ratio for this system. Al III was 
detected only in the components at v= —197 and v= —162 km 
s~^, and not in the components where other ions were the bulk of 
the absorption from other ions was seen, such as the v=2, v=38, or 
v=74 km s^^ components. We did not detect any Zn II lines in this 
system at S/N'~23 in the region, and place an upper limit on the Zn 
metallicity of [Zn/H]< —1.14. The Fe II abundance is similarly 
low for this system, with [Fe/H]= — 1.33±0.07. We show velocity 
plots of several lines in Figure 5. The results of the profile fitting 
analysis are found in table 8, with the rest frame equivalent widths 
for this system in table 2. 



4.3 Q1215-0034 (z<,m=2.691) 



(System A:za6s=1.5543): This fairly hi gh redshift QSO harbors 
a sub -DLA with log Nh i =19.56±0.02 jRao, Tumshek. fc Nesto3 
I2OO6D . Fourteen components were used in the profile fit for this 
system, spanning -^390 km s"'^ in velocity space. Several of the 
components were detected only in the Mg II AA 2796, 2803 lines. 
We detected lines from Mg II AA 2796, 2803 and Fe II AA 2344, 
2374, 2382, 2586, and 2600. The Al II A 1670 and Al III AA 1854, 
1862 lines were covered, but are blended with features inside the 
Lyman-a forest. The Mg I A 2852 line was blended with telluric 
features at ~7284 A. We did not detect Zn II at S/N ~30 in the 
region, but we could place an upper limit of [Zn/H]< —0.56, and 
based on the Fe II lines an abundance of [Fe/H]= — 0.64±0.02. 
Velocity plots of lines of interest are shown in Figure 4. The results 
of the profile fitting analysis are given in table 7, with the rest frame 
equivalent widths for the lines in this system in table 2. 



4.4 Q1220-0040 (Zem=1.411) 



sub-DLA with log 
spectrum of this QSO 



(System A:zai,s=0.9746): There is a 
Nh I =20.20±0.07 in the UV 

( iRao, Tumshek. fc Nesto3 I2OO6I) . We detect lines of Al III 
Mg I, Mg II, and Fe II in this sub-DLA. A total of 9 components 



4.5 Q1228-1-1018 (z^™=2.305) 

(System A:Zabs=0.9376): This moderately high redshift QSO 
has a sub-DLA in its s pectrum with log Nh i =19.41 ±0.02 
dRao, Tumshek. fc Nestoill2006.) . We detected strong lines of Mg 
I, Mg II, and Fe II in this sub-DLA. A total of 9 componenets 
were used in the Voigt profile fits. Due to the high redshift of the 
background QSO, several lines were blended with Lyman-a forest 
clouds, such as Al III A A 1854, 1862 and Si II A 1808. The Zn II 
AA 2026, 2062 lines also fell inside the Lyman-a forest, but lay be- 
tween strong absorption features, especially the Zn II A 2026 line. 
There is a strong forest line in the blue end of the expected position 
of the Zn II feature, as can be seen in the velocity plots in Figure 6, 
however the components which are most likely to hold Zn II are at 
v=2, v=26, and v=51 km s~^, which are unblended. Nonetheless, 
no Zn II is seen, and we place an upper limit on the Zn metallicity 
of [Zn/H]< —0.37 with S/N~26 in the region. The Fe abundance 
for this sub-DLA is high compared to typical DLA Fe abundances, 
with [Fe/H]= -0.30±0.02. The results of the profile fitting anal- 
ysis are given in table 9, with the rest frame equivalent widths for 
the lines in this system given in table 2. 



Table 2. Rest-frame equivalent widths of key metal lines. Values and la errors are in units of mA. 



QSO 


^abs 


Mgl 


Mgll 


Mgll 


Al 11 


Aim 


Aim 


Sill 


Sill 


Call 


Ca 11 


Cr 11 






2852 


2796 


2803 


1670 


1854 


1862 


1526 


1808 


3933 


3969 


2056 


Q 1037+0028 


1.4244 


a 


2443±19 


2054±32 


677±15 


169±12 


87± 8 


746±49 


60±9 


- 


- 


<7 


Q1054-0020 


0.8301 


350±15 


1180±14 


1077±12 


- 


415±18 


287±28 


- 


- 


- 


- 


<11 




0.9514 


131±22 


769±19 


542±27 


- 


<22 


<11 


- 


<11 


- 


- 


<6 


Q1215-0034 


1.5543 


- 


1987±21 


1543±26 


- 


- 


- 


- 


- 


- 


- 


<10 


Q 1220-0040 


0.9746 


152±20 


1919±32 


1688±32 


- 


71±22 


40±22 


- 


- 


- 


- 


<9 


Q1228+1018 


0.9376 


245±38 


1678±40 


1457±23 


- 


- 


- 


- 


- 


- 


- 


<6 


Q1330-2056 


0.8526 


293±66 


2413±64 


1807±41 


- 


75±24 


26± 19 


- 


- 


- 


- 


<6 


Q1436-0051 


0.7377 


583±14 


1113±10 


1081±13 


- 


- 


- 


- 


- 


396±15 


- 


<20 




0.9281 


233±21 


1105±40 


867±33 


- 


220±22 


156±22 


- 


46±9 


122±17 


- 


<16 


Q1455-0045 


1.0929 


363±21 


1659±10 


1586±11 


804±14 


515±15 


281±17 


- 


- 




- 


21 ±7 


QSO 


^abs 


Mn 11 


Mn 11 


Mn 11 


Fell 


Fe 11 


Fe 11 


Fell 


Fell 


Fell 


Znll 


Znll 






2576 


2594 


2606 


2260 


2344 


2374 


2382 


2586 


2600 


2026 


2062 


Q 1037+0028 


1.4244 


75±20 


82±18 


40± 13 


93±21 


963±61 


507±21 


1442±55 


839±51 


1427±41 


<21 


<15 


Q1054-0020 


0.8301 


42±11 


39±8 


<8 


<7 


586±7 


277±9 


805±6 


548±6 


831 ±6 


<10 


<15 


Q 1054-0020 


0.9514 


<18 


<10 


<13 


<6 


198±15 


45±9 


420±18 


208±27 


453±27 


<9 


<13 


Q1215-0034 


1.5543 


<9 


<8 


<7 


24±5 


507±14 


267±9 


861 ±20 


449±10 


896±24 


<8 


<10 


Q 1220-0040 


0.9746 


<19 


<22 


<17 


<8 


624±10 


278±10 


1049±32 


566±25 


932±22 


<9 


<9 


Q1228+1018 


0.9376 


<29 


<20 


<18 


64±12 


692±12 


417±22 


989±9 


660±19 


1032±26 


<9 


<8 


Q1330-2056 


0.8526 


<9 


<8 


<8 


<6 


271±51 


58± 43 


803±43 


204±41 


784±40 


<16 


<6 


Q1436-0051 


0.7377 


207±17 


134±14 


85±15 


116±24 


757±17 


580±16 


847±9 


785±11 


926±16 


108±12 


101 ± 35 




0.9281 


43±16 


<22 


102±32 


25±7 


341±16 


189±13 


501±13 


307±36 


471 ±23 


26±7 


<21 


Q1455-0045 


1.0929 


29±6 


50±14 


<19 


<12 


902±11 


470±22 


1135±14 


874±24 


1226±21 


<15 


<8 



"Blank entries indicate that the line was not measured due to lack of wavelength coverage, extremely low S/N at the wavelength extremes, or blends with 
telluric features. 



Table 3. Total column densities for tiie systems observed. Values are logarithmic, in units of cm 





^abs 


log Nh I 


Mgl 


Mg n 


A 1 TT 
Al 11 


A 1 TTT 
Al 111 


TT TT 

5>i 11 Ca 11 


Ti n 


TT 

Cr 11 


Mn n 


Fe n 


Zn n 


Q1037+0028 


1.4244 


20.04±0.12 


_ 


> 15.48 


>14.15 


13.19±0.02 


15.09±0.03 


_ 


< 12.26 


12.53±0.05 


14.84±0.02 


<12.04 


AOD 








> 14.44 


> 13.60 


13.15±0.03 


15.05±0.05 






12.57±0.11 


14.96±0.09 




Q1054-0020A 


0.8301 


18.95±0.18 


12.55±0.01 


> 14.83 




13.59±0.06 






< 12.47 


12.27±0.07 


14.33±0.01 


<11.76 


AOD 






12.54±0.02 


> 14.30 




13.65±0.04 








12.31±0.10 


14.35±0.01 




Q1054-0020B 


0.9514 


19.28±0.02 


12.04±0.03 


> 13.72 




<12.11 


< 14.25 


< 12.36 


< 12.20 


< 11.92 


13.66±0.01 


< 11.70 


AOD 






12.04±0.07 


> 13.59 














13.49±0.08 




Q1215-0034 


1.5543 


19.56±0.02 




>15.34 










< 12.37 


<11.63 


14.39±0.01 


<11.63 


AOD 








> 14.25 














14.35±0.11 




Q 1220-0040 


0.9746 


20.20±0.07 


12.12±0.03 


>14.81 




12.62±0.06 




< 12.47 


< 12.38 


<11.95 


14.34±0.02 


<11.69 


AOD 






12.10±0.05 


> 14.37 




12.68±0.12 










14.35±0.01 




Q1228+1018 


0.9376 


19.41±0.02 


12.32±0.04 


> 14.61 








<11.65 


<12.21 


<12.14 


14.58±0.01 


<11.67 


AOD 






12.32±0.05 


> 14.36 














14.55±0.02 




Q1330-2056 


0.8526 


19.40±0.02 


12.33±0.03 


> 14.21 




12.59±0.05 




< 11.45 


<12.18 


<11.61 


13.80±0.01 


<11.96 


AOD 






12.37±0.09 


>14.12 




12.52±0.15 










13.76±0.07 




Q1436-0051A 


0.7377 


20.08±0.11 


> 12.92 


> 14.42 






12.83±0.02 


<11.71 


<12.71 


13.01±0.02 


14.94±0.02 


12.67±0.05 


AOD 






> 12.90 


> 14.26 






12.79±0.01 






13.05±0.03 


14.92±0.02 


12.80±0.05 


Q1436-0051B 


0.9281 


<18.8 


12.49±0.04 


>15.12 




13.40±0.03 


15.02±0.05 12.28±0.04 


< 12.71 


< 12.59 


12.45±0.05 


14.20±0.01 


12.29±0.06 


AOD 






12.45±0.03 


>14.11 




13.31±0.03 


14.97±0.07 12.15±0.05 






12.32±0.07 


I4.25±0.02 


12.23±0.05 


Q1455-0045 


1.0929 


20.08±0.06 


12.53±0.02 


> 14.92 


>13.83 


13.65±0.02 


14.64±0.10 


< 12.44 


12.40±0.13 


12.09±0.15 


14.57±0.01 


<11.91 


AOD 






12.51±0.02 


> 14.45 


>13.71 


13.67±0.01 


14.84±0.12 




12.75±0.17 


12.20±0.10 


14.61±0.02 





Sub-DLAs and strong Lyman-limit systems at 1.5. 1 



Table 4. Column densities for the z =1.4244 absorber in Q1037+0028 with Nh i =20.04±0.12. Velocities and be// 
values are in units of km s~^, and column densities are in units of cm~^. Errors in the column densities are the la 
values. 



V 




Mgn 


Ain 


Aim 


Sin 


Mnn 


Fen 


-158 


6.9 


> 3.14E13 


> 2.34E12 


(6.59±1.40)E11 


(4.26±0.59)E13 




(1.63±0.13)E13 


-118 


5.9 


> 8.37E13 


> 2.16E12 




(5.69±0.98)E13 




(3.04±0.21)E13 


-91 


6.7 


> 3.91E12 


> 4.06E11 




(6.52±0.20)E12 




(2.92±0.76)E12 


-57 


13.5 


> 1.65E13 


> 1.28E12 




(3.65±0.46)E13 




(1.65±0.14)E13 


-40 


6.9 


> 9.77E11 


> 1.47E12 




(2.98±0.60)E13 




(1.69±0.17)E13 


-19 


10.0 


> 6.90E13 


> 4.60E12 


(1.37±0.17)E12 


(7.21±0.93)E13 




(3.92±0.21)E13 


11 


8.1 


> 1.23E14 


> 5.88E13 


(1.12±0.05)E13 


(7.10±0.61)E14 


(2.07±0.29)E12 


(3.08±0.25)E14 


44 


7.9 


> 2.69E15 


> 6.80E13 


(2.42±0.18)E12 


(2.53±0.48)E14 


(1.33±0.26)E12 


(2.36±0.17)E14 


100 


11.2 


> 7.16E12 


(4.64±0.98)E11 




(6.79±2.23)E12 




(5.65±0.89)E12 


147 


12.6 


> 1.56E13 


(9.89±1.16)E11 




(2.43±0.32)E13 




(1.33±0.11)E13 


183 


4.6 


> 4.57E12 













Table 5. Same as table 3, but for the Zai,s=0.8301 system in Q1054-0020 with Nh i =18.95±0.18. 



V 


be// 


Mgl 


Mgn 


Aim 


Mnn 


Fen 


-25 


9.1 


(4.53±0.41)E11 


>1.07E14 


(9.41±1.86)E12 




(3.07±0.08)E13 


2 


10.3 


(1.53±0.08)E12 


>2.00E14 


(1.80±0.43)E13 


(8.30±1.68)E11 


(7.96±0.20)E13 


19 


8.1 


(6.91±0.58)E11 


> 8.82E12 




(4.12±1.50)E11 


(2.63±0.10)E13 


49 


13.4 


(8.89±0.51)E11 


> 3.59E14 


(1.15±0.18)E13 


(6.09±1.73)E11 


(7.94it0.15)E13 


104 


15.8 




(1.23±0.12)E12 









Table 6. Same as table 3, but for the Zai,s=0.9514 system in 
Q1054-0020 witii Nh i =19.28±0.02. 



V 


be// 


Mgl 


Mgn 


Fen 


-447 


8.0 




(1.15±0.11)E12 




-90 


5.1 


(1.71±0.28)E11 


>6.74E12 


(4.86±0.27)E12 


-71 


4.6 




>1.84E12 




-13 


10.0 


(2.23±0.31)E11 


>1.24E13 


(9.84±0.34)E12 


49 


3.5 


(7.26±2.45)E10 


>1.43E12 


(1.71±0.23)E12 


72 


8.0 


(5.24±0.41)E11 


>2.64E13 


(2.66±0.06)E13 


87 


10.3 


(1.07±0.33)E11 


>3.07E12 


(2.69±0.33)E12 



Table 7. Same as table 3, but for the Zoi,s=1.5543 
system m Q1215-0034 witii Nh i =19.56±0.02. 



V 


be// 


Mgn 


Fen 


-271 


8.0 


(1.42±0.09)E12 




-217 


6.7 


(1.84±0.10)E12 




-188 


9.7 


>7.931E12 


(2.52±0.41)E12 


-154 


9.0 


>1.10E13 


(2.28±0.40)E12 


-131 


7.4 


>1.04E12 




-103 


13.9 


>1.77E12 




-75 


8.7 


>5.48E12 


(5.22±0.44)E12 


-45 


10.0 


>3.30E13 


(1.75it0.06)E13 


-3 


9.7 


>2.06E15 


(9.13±0.36)E13 


14 


15.5 


>6.78E13 


(1.18±0.03)E14 


42 


10.9 


>1.57E13 


(6.80±0.61)E12 


70 


10.9 


(2.92±0.13)E12 




102 


5.4 


(1.16±0.09)E12 




120 


6.2 


(1.87±0.10)E12 


(1.77±0.37)E12 
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Table 8. Same as table 3, but for the z„5s=0-9746 system in Q1220- 
0040 with Nh I =20.20±0.07. 



V 




Aim 


Mgl 


Mgll 


Fell 


-197 


4.9 


(1.50±0.33)E12 




>6.51E12 


(1.03±0.27)E12 


-162 


5.6 


(2.70±0.42)E12 




>1.53E13 




-135 


11.1 






>4.45E12 




-100 


4.6 






>4.19E12 


(2.08±0.29)E12 


-45 


8.0 






>3.10E12 


(1.52±0.32)E12 


2 


14.6 




(6.71±0.52)E11 


>4.95E14 


(1.22±0.02)E14 


38 


15.8 




(3.91±0.48)E11 


>5.58E13 


(4.91±0.11)E13 


75 


15.2 




(2.59±0.44)E11 


>5.86E13 


(4.08±0.09)E13 


106 


9.0 






>6.07E12 


(2.60±0.36)E12 



Table 9. Same as table 3, but for the Zat,s =0.9376 system in 
Q1228+1018 with Nh i =19.41±0.02. 



V 




Mgl 


Mgll 


Fell 


-44 


20.2 




(3.05±0.27)E12 


(3.05±0.58)E12 


-15 


8.9 


(2.16±0.69)E11 


>6.96E13 


(1.86±0.08)E13 


2 


9.2 


(1.18±0.66)E11 


>2.46E13 


(3.97±0.14)E13 


26 


9.9 


(9.72±1.02)E11 


>1.64E14 


(2.29±0.05)E14 


51 


11.1 


(4.93±0.78)E11 


>6.42E13 


(3.95±0.12)E13 


77 


8.0 


(1.88±0.61)E11 


>5.86E13 


(4.31±0.13)E12 


109 


7.4 




>1.39E13 


(6.54±0.45)E12 


126 


4.9 


(1.04±0.36)E11 


>8.37E12 


(1.72±0.38)E12 


155 


14.9 




(4.17±0.25)E12 


(1.92±0.48)E12 



4.6 Q1330-2056 (Zem=1.165) 

(System A:za6s=0.8526): In the spectrum of this relatively bright 
(mv = 17.0) QSO, there is a sub-DLA at Zabs= .8526 with log 
Nh I =19.40±0.02 jRao. Turnshek. fc Nestorll2006l) . We detected 
strong absorption from Al III, Mg I, Mg II, and Fe II. The Al II A 
1670 line was below the covered wavelengths, so the Al III/Al II 
ratio could not be determined for this system. This system required 
15 components in the Voigt profile models for an adequate fit. The 
components spanned ~650 km s~^ in velocity space. The Zn II 
A 2026 line fell near a bad column in the CCD on the blue side. 
However, the components where Zn II is most likely, such as the 
components with the strongest Mg I absorption, are outside of the 
affected area. We have highlighted the affected area in grey in the 
plot of the Zn II A 2026 line in Figure 7. No Zn II was detected, 
and an upper limit of [Zn/H]< —0.07 was placed on the system at 
S/N~18 in the region. The Fe abundance was low, with [Fe/H]< 
— 1.07±0.02. We give velocity plots of several lines in Figure 7. 
The results of the profile fitting analysis are shown in table 10, with 
the rest frame equivalent widths for the lines in this system given 
in table 2. 

4.7 Q1436-0051 (Zem=1.275) 

(System A:Zni,s=0-7377): This is a s ub-DLA system with log Nh i 
=20.08±0.11 ( Rao, Turns hek. fc Ne stor 2006). We detect absorp- 
tion Hnes from Mg I, Mg II, Ca II, Mn II, Fe II, and Zn II in 
system A. A total of 8 components were used to fit the absorp- 
tion profiles. There is also a C IV system at Zn6ji = 1.2699, with the 
C IV A 1550 line close to the expected position of the Zn II A 



2026 line in the sub-DLA at Za(,s =0.7377, as can be seen in the 
velocity plots shown in Figure 8. The C IV AA 1548, 1550 lines 
were fit simultaneously, and later held fixed while fitting the Zn 
II AA 2026, 2062 lines simultaneously. The C IV system shows a 
simple structure with only 2 components needed to fit the profile, 
with no absorption in the region where the majority of the Zn II A 
2026 absorption is. In Figure 8, the C IV contribution to this line 
is highlighted in grey. The Mg I contribution to the Zn II A 2026 
line, based on the fit to the Mg I 2852 line, was also included and 
held fixed. This was the only case of an appreciable contribution 
of Mg I to the blended Zn II -I- Mg I A 2026 line that was seen. 
This system has a high metallicity with [Zn/H]= — 0.05±0. 11, and 
[Fe/H]= -0.61±0.1 1. We also detect strong Mn II A A 2576, 2594, 
2606 lines with [Mn/H]= -0.57±0.11 and [Mn/Fe]=-^0.04±0.03. 
The Ca II A 3934 line was detected with Wo=396±15 mA, but the 
Ca II A 3969 line was blended with telluric features. The Si II A 
1808 line was not observable for this sub-DLA, so no information 
on aZp/ia-enhancement is available for this system. We show the 
results of the profile fitting analysis in table 1 1 , and give the rest 
frame equivalent widths for this system in table 2. 

(System B:Zabs=0-9281): This is a Lyman-limit system with 
log Nh I <18.8 (Rao, Turnshek, & Nestor 2006). Several lines are 
detected in this system, such as Mg I, Mg II, Al III, Si II, Ca II, 
Mn II, and Fe II. A total of 6 components were used in the Voigt 
profile fits for this system, although Si II, Ca II, and Zn II were 
only detected in the v ~9 km s~^ component. This system has a 
large amount of Al III and likely has a high degree of ionisation, 
although the Al II A 1670 line was not in the wavelengths covered. 
The Ca II A 3969 line was again blended with telluric features. The 
Fe II A 2344 line is partly contaminated by the Fe II A 2600 line 
from system A. 
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Figures 1 and 2: Velocity plots for several lines of interest in the z — 1 .4244 system (top) in the spectrum of Q 1037+0028 and the z^bs 

=0.8301 

system in Q1054-0020 (system A, bottom). The solid green line indicates the theoretical profile fit to the spectrum, and the dashed red line is 
the continuum level. The vertical dotted lines indicate the positions of the components that were used in the fit. In the cases of the Zn 11 AA 
2026,2062 lines, which have other lines nearby, the long dashed vertical lines indicate the positions of the components for Mg I (former case), 
and Cr II (latter case). 
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Figures 3 and 4: Same as Fig lire 1, but for the z^^fj^— 0.9514 system in Q1054-0020 (system B, top), and the z^^g— 1.5543 system in Q1215- 
0034 (bottom). For the ^abs~l-^543 system in Q1215-0034, the area shaded in grey in the panel displaying the Zn II A 2026 line is an 
unidentified feature, not due to absoiption from Zn II. This may be due to Mg I 2026, however, the Mg I 2852 line was blended with telluric 
features, and we could therefore not estimate the Mg I contribution to this line. 
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Figures 5 and 6: Same as Figure 1, but for the Zabs=0.9746 system in Q1220-0040 (top), and the Zabs=^-9316 system in Q1228+1018 
(bottom). The ai'eas shaded in grey in the Zn II AA 2026, 2062 panels for the 7^^^,3=0.9376 system in Q1228+1018 represent absorption from 
lines inside the Lyman-a forest, and are not due to absorption from Zn II. There is a small, unidentified feature near the Fe II A 2374 line, 
which has also been shaded in grey. 
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Figures 7 and 8: Same as Fig ure 1, but for the z^jj^— 0.8526 system in Q1330-2056 (top), and the z^^g— 0.7377 system in Q1436-0051 (system 
A, bottom). For the z„5s=0.8526 system in Q1330-2056, in the region around the expected position of Zn II A 2026, there is a bad column in 
the blue CCD that caused the feature seen in the figure. The part of the profile that is affected by the bad column is shaded in grey to highlight 
this defect. For the Za,,3=0.7377 system in Q1436-0051, we have highlighted the contribution of the C IV A 1550 line at Zni,3=1.2699 to the 
Zn II 2026 profile in grey. See § 4.7 for a discussion of this blend. 
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Figures 9 and 10: Same as Fig ure 1, but for the ^.ahs 

=0.9281 system in Q1436-0051 (system B, top), and the 7^1,3=1.0929 system in Q1455- 

0045 (bottom). 
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Table 10. Same as table 3, but for the z^bs=0.8526 system in Q1330-2056 with 
Nh I =19.40±0.02. 



V 




Al III 


Mgl 


Mg II 


Fe II 


-439 


6.2 


(6.39±1.74)E11 


- 


>9.07E12 


(2.23±0.37)E12 


-365 


9.9 


(1.91±0.23)E12 


(1.80±0.40)E11 


>3.09E13 


(8.13±0.50)E12 


-295 


4.3 






(1.44±0.28)E12 




-262 


6.4 




(1.03±0.35)E11 


(2.23±0.30)E12 




-192 


9.7 


(6.60±1.96)E11 


(2.84±4.28)E11 


>2.05E13 


(6.50±0.47)E12 


-121 


7.2 






(2.28±0.33)E12 




-101 


10.4 




(2.07±0.43)E11 


>5.25E12 


(1.26±0.42)E12 


-74 


12.8 






>3.49E12 


(1.56±0.46)E12 


-10 


10.9 






(2.31±0.31)E12 




21 


9.2 




(2.97±0.43)E11 


>2.91E13 


(9.65±0.52)E12 


46 


10.4 




(1.17±0.41)E11 


>3.83E12 




67 


11.4 




(1.09±0.41)E11 


>6.20E12 


(3.02±0.47)E12 


99 


7.5 




(4.30±0.48)E11 


>2.74E13 


(1.75±0.07)E13 


119 


7.9 




(1.58±0.41)E11 


>7.84E12 


(7.74±0.50)E12 


136 


6.4 


(6.36±1.78)E11 


(1.12±0.38)E11 


>1.04E13 


(4.31±0.42)E12 



This system appears to have a super-solar abundance, with 
[Zn/H]>+0.86 and based on Fe II, the metallicity is [Fe/H]> 
—0.07. As was mentioned in § 4.2, the [Fe/H] may be lower 
than this due to ionisation corrections. This system also shows a 
super-solar abundance in both Mn with [Mn/H]> -1-0.15 and Si with 
[Si/H]>-l-0.68. If the Mg I contribution to the blended Zn II -I- Mg 
1 A 2026 line is substantially more than what we have estimated 
based on the Mg I A 2852 line, then the Zn II metallicity could be 
overestimated. However, the component structure of Mg I A 2852 
does not suggest that an Mg I A 2026 blend could affect the de- 
tected dominant Zn II A 2026 component. With such a low H I 
column density, it is possible that a substantial ionisation correc- 
tion factor is needed for this system. The Al III/Al II ratio which 
can be used to estimate ionisation, could not be determined for this 
system as the Al II A 1670 was not covered. Even with —1.0 dex 
ionisation corrections lowering the abundances, this system would 
still be well above the average metallicity of DLAs. Figure 9 gives 
velocity plots of several lines of interest. We give the results of the 
profile fitting analysis in table 12, along with the rest frame equiv- 
alent widths for the lines in this system in table 2. 



4.8 Q1455-0045 (Zem=1.378) 

(System A:zabs=1.0929): There is a sub-DLA system 
with log Nh i =20.08±0.06 in the spectrum of this QSO 
( iRao, Tumshek. &: Nestoj [2006h . We detected strong absorption 
lines from Mg I, Mg II, Al II, Al II, Si II, and Fe II in this 
system. A total of 10 components were used in the Voigt profile 
fits. No Zn II AA 2026, 2062 lines were detected at S/N~27 
in the region. We could constrain the metallicity based on Zn 
as [Zn/H]< —0.80. The abundances based on the Fe II and Si 
n lines are [Fe/H]= -0.98±0.06 and [Si/H]= -0.98±0.12. 
The Mn II A 2576 was also detected with [Mn/H]=-1.43 and 
[Mn/Fe]= -0.51±0.15 We did detect both the Al II and Al III 
lines, with the ratio of column densities Al III/Al II < —0.18. Weak 
features of Cr II A 2056 were also detected with Wo=21±7 mA. 
Only one component at ~64 km s~^ was above a 3a significance 
level with [Cr/Fe]= — 0.35±0. 13, although another component at 



~ 44 km s~^ was at the ~2(t level. The Cr II column density was 
found to be higher via the AOD method, with [Cr/Fe]=0.0. Figure 
10 shows velocity plots of several lines of interest. We show the 
results of the profile fitting analysis in table 13, and the rest frame 
equivalent widths for this system in table 2. 



5 DISCUSSION AND CONCLUSIONS 

The abundances for the observed systems are given in table 14. We 
have used the total column densities (i.e. the sum of the column 
densities in the individual components of a system that were de- 
termined via profile fitting method) along with the total Nh i as 
given in table 1, to determine the abundances of these systems. 
It would be interesting to determine the abundances in the indi- 
vidual components, as they may possibly vary from component to 
component. This however is impossible at the present time due to 
the lack of high resolution UV spectrographs that would be neces- 
sary to determine Nh i in individual components. We have ignored 
any ionisation corrections while determining these abundances, and 
have assumed the first ions to be the dominant ionisation species of 
the elements for which these abundances have been determined, 
namely Zn, Fe, Mn, Cr, and Si. We discuss th e issue of ionisation 
further below. Solar systems abundances from lLodderj ( l2003h are 
also given in table 14. The abundances for Ca II are likely only 
lower limits, as the ionisation potential of 11.868 eV for Ca"*" is 
lower than ionisation potential of H" (13.6 eV). 

Relative abundances of various elements are also given in ta- 
ble 14, with the column densities determined from the profile fitting 
analysis. We give the metallicities of these systems, along with the 
[Zn/Fe] ratio, which is often used as an indicator of dust deple- 
tion. We also provide the ratios of [Si/Fe], [Ca/Fe], [Cr/Fe], and 
[Mn/Fe]. Finally, we provide ratios of the column densities of the 
adjacent ions Al III/Al II and Mg II/Mg I as well as Mg II/ Al III 
and Fe II/ Al III, any of which may provide information about the 
ionisation in these systems. By matching the observed ratios of the 
ions mentioned above to those generated through grids of photoion- 
isation models, one can constrain the levels of ionisation, and any 
ionisation correction factors for the observed abundances. 

Based on grids of Cloudy models, is has previously been 
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Table 11. Same as table 3, but for the z„i,s=0.7377 system in Q1436-0051 with Nh i =20.08±0.11. 



15 



V 




Mgl 


Mg 11 


Call 


Mnll 


Fell 


Znll 


-49 


3.6 




>1.00E13 




(4.40±1.()7)E11 






-34 


13.2 


(6.65±0.76)E11 


>8.68E13 


(4.18±1.37)E11 


(2.18±0.21)E12 


(1.56±0.06)E14 




-25 


3.5 


(2.67±0.81)E11 






(6.66±1.53)E11 


(2.09±0.12)E13 




-11 


11.1 


>2.51E12 


>5.74E13 


(2.69±0.17)E12 


(4.11±0.19)E12 


(3.82±0.33)E14 


(2.17±0.49)E12 


19 


11.5 


>3.48E12 


>7.13E13 


(2.61±1.68)E12 


(2.47±0.16)E12 


(2.70±0.16)E14 


(1.74±0.17)E12 


35 


5.6 


>1.31E12 


>2.26E13 


(9.22±1.25)E12 


(4.29±I.12)E11 


(3.00±0.34)E13 


(7.16±1.36)E11 


47 


5.7 


(1.15±0.36)E11 


>1.23E13 






(6.94±0.71)E12 




62 


8.0 




(4.03±0.21)E12 






(2.24±0.52)E12 





Table 12. Same as table 3, but for the z„6s=0.9281 system in Q 1436-0051 withNni <18.8. 



Mgl 



Mg 11 



Al 111 



Sill 



Ca 11 



Mn 11 



Fell 



Znll 



-70 6.8 

-44 6.6 

-14 2.7 

9 10.7 (2.79±0.28)E12 

36 7.3 

60 8.0 (2.86±0.75)E11 



(1.70±0.37)E12 
(3.00±0.51)E12 
>1.09E15 
>9.99E13 
>1.61E13 
>1.15E14 



(4.62±0.93)E12 
(1.83dz0.14)E13 
(1.03±0.33)E12 
(1.17±0.33)E12 



(1.04±0.13)E15 (1.89±0.I6)E12 



(4.64±1.82)E11 
(2.36±0.27)E12 



(8.56±1.00)E12 
(1.30±0.05)E14 
(6.91±0.70)E12 
(1.26±0.08)E13 



(1.97±0.27)E12 



shown that there appears to be little need for ionisation coiTec- 
tions in most sub-D LA systems jPessauges-Zavadskv et al. I l2003l : 
iMeiring et al.ll2007l) . We have therefore not applied any ionisation 
corrections to to the abundances in table 14. We do note however, 
that the z„bs=0.8301 system with Nh i =18.95±0.18 in Q1054- 
0020 and the z^bs =0.9281 system with Nh i <18.8 in Q 1436-0051 
both show possible signs of high levels of ionization, with Fe II/ Al 
111=0.74 and Fe II/ Al 111=0.80 respectively. Substantial ionisation 
corrections could be necessary. Cloudy models for the Zai,s=0.9281 
system in Q 1436-005 1 indicate a positive ionization correction fac- 
tor for Zn in this system of ~-l-0.2 to -1-0.3 dex for —6 <log U< —3 
(indicating that [Zn II/H I] underestimates the true abundance by 
this amount). We caution that these results may be inaccurate due 
to uncertainties in the atomic data for Zn. The ionisation correction 
factors for Si, Mn, and Fe are all negative over this range of ionisa- 
tion parameters with values of ~ —0.3 to —0.5 dex for Mn and Fe, 
and values ranging from '--^O.O to ~ —1.25 for Si. Based on these 
models, it appears that [Fe/Zn] is suppressed as one goes to lower 
and lower Nh i values, the true depletion is actually much higher 
than it appears based on Zn II and Fe II. Similar values were derived 
for the Zabs=0.830I system with Nh i =I8.95±0. 18 in Q1054-0020 
based on Cloudy models. As no adjacent ions that could constrain 
the ionisation parameter are available for these systems, we leave 
the abundances unaltered in table 14. We prefer to leave the correc- 
tions to a later time when more constraints are available (i.e. when 
we have UV spectra with coverage of Fe III, S III, Si IV, etc and 
better atomic data become available). 

Variations by component of the relative abundances of ele- 
ments can also provide informati on about the local star forma- 
tion history and possibly ionization. IProchaska et alj ( l2002h noticed 
a remarkably small (~0.3 dex) component-by-component disper- 
sion in [Zn/Fe] over a wide range of metallicities. As we only 
have Zn detections for two of the systems here, we instead look 
at [Mn/Fe], which is likely less affected by dust depletion (Mn and 
Fe have similar condensation temperature) and more affected by 
nucleosynthetic enrichment. Unfortunately, other ions such as Mg 



II and Al II are too saturated to determine accurate column den- 
sities in most of the individual components. The Za6s=1.4244 sys- 
tem in QI037-I-0028 and the z„i,s=0.830I system in Q1054-0020 
show similarly low dispersions of [Mn/Fe], but interestingly the 
high metallicity Z£i(,s=0.7377 and Za6s=0.928I systems of QI436- 
005 1 show higher spreads of [Mn/Fe] of ~0.6 and ■-^0.5 dex respec- 
tively. The metal rich sub-DLA observed in lPeroux et alj ( l2006al) 
also showed a large dispersion in [Mn/Fe]. The z^bs =0.7377 sys- 
tem in Q 1436-005 1 also shows a 0.6 dex range of [Zn/Fe] in the 
components. One possible explanation for these larger dispersions 
could be enrichment by recent supemovae explosions in which the 
gas has not had time to sufficiently mix. 

Zn II was detected in only 2 of the systems, and upper lim- 
its on the Zn abundance were placed on the rest of the systems. 
The Cr II A 2056 line was detected only in the Za6s=1.0929 sys- 
tem in Q 1455-0045. Due to the weakness of the Lyman-a line 
in the Za53=0.928I system in QI436-005I, only an upper limit 
could be placed on the neutral hydrogen column density, and there- 
fore lower limits on the abundances in this system. This sys- 
tem shows an unusual abundance pattern, with high metallicity 
([Zn/H]>-l-0.86, [Fe/H]> -0.07, and [Si/H]>-l-0.68) and high de- 
pletion ([Zn/Fe] =-1-0. 94). This system also showed an overabun- 
dance of Mn, with [Mn/Fe] =-1-0.22. As Mn is an odd element, 
and Fe is even, [Mn/Fe] might be below solar in material whose 
source is dominated by the products of Type II supernovae, but 
likely not higher if differential dust depletion effects are negligi- 
ble. In the other case of a Zn detection, the system at Za(,s=0.7377 
in QI436-005I, Mn is also slightly more abundant r elative to Fe, 
with [Mn/Fe]=-l-0.04. lLedoux. Bergeron, fc PetitieanI 12002) found 
a large dispersion in [Mn/Fe] based on a sample of DLA absorbers 
and concluded that this large range of values not seen in other abun- 
dance ratios such as [Cr/Fe] or [Si/Fe] must be due to a combination 
of both depletion and nucleosynthetic effects. 

A similar over-abu ndance of Mn relative to Fe was also seen 
by IPeroux et alj i2006ah in another sub-DLA at Za(,s=0.716 with 
[Zn/H]=-l-0.61 and [Mn/Fe]=-l-0.I4. This could perhaps be due to 
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Table 13. Same as table 3, but for the z„i,s=1 0929 system in the spectrum of Q 1455-0045 with Nh i =20.08±0.06. 



V 




Mgl 


Mg II 


AlII 


Al III 


Sin 


Cr II 


Mn II 


Fell 


-95 


8.8 




(1.00±0.09)E12 




(4.92±1.51)E11 










-60 


9.8 


(3.01±0.50)E11 


>2.02E14 


>7.21E12 


(2.87±0.22)E12 








(2.15±0.09)E13 


-42 


7.5 


(7.42±0.78)E11 


>8.96E12 


>3.09E12 


(2.46±0.23)E12 






(1.23±0.40)E12 


(3.64±0.20)E13 


-30 


8.2 


(4.12±0.59)E11 


>1.33E14 


>7.51E12 


(2.37±0.22)E12 








(5.02±0.22)E13 


3 


10.3 


(1.81±0.55)E11 


>1.10E14 


>1.62E13 


(1.06±0.05)E13 








(4.73±0.19)E13 


15 


7.1 


(2.95±0.65)E11 


>1.05E13 


>5.51E12 


(4.69±0.50)E12 








(2.24±0.20)E13 


26 


8.1 




>2.20E13 


>1.80el2 


(2.03±0.30)E12 








(1.51±0.13)E13 


45 


11.1 


(4.89±0.59)E11 


>7.02E13 


>1.67E12 


(9.50±0.44)E12 








(6.78±0.26)E13 


64 


7.9 


(6.30±0.71)E11 


>5.66E13 


>3.05E12 


(7.52±0.46)E12 




(2.51±0.71)E12 




(6.66±0.38)E13 


77 


8.8 


(2.17±0.53)E11 


>2.18E14 


>6.12E12 


(2.56±0.25)E12 


(4.40±0.10)E14 






(4.65±0.21)E13 



heavy depletion of Fe compared to Mn onto dust grains, unmixed 
gas that has recently been expelled through stellar winds and su- 
pemovae, or an intrinsically non-solar abundance pattern. Cloudy 
models for this system indicate that the corrections for [Mn/Fe] 
due to ionisation are small (< 0.10 dex for —6 <log U< —1). 
In the other 3 cases where we have detections of Mn and Fe, 
but no Zn, Mn is slightly under- abundant relative to Fe, with 
<[Mn/Fe]>=— 0.31. The systems in our sample have an average 
Fe abundance of <[Fe/H]>=— 0.68, which is higher than for typi- 
cal DLA systems (<[Fe/H]>DLA~ —1.5). This may be due to less 
dust depletion or truly higher abundances in these systems, or to 
undue influence of ionisation effects in the low Nh i systems. 

To determine an average [Zn/H] abundance of the absorbers in 
this sample, the spectra were shifted to rest wavelengths and aver- 
aged. Any obvious non-associated features in the region around the 
Zn II A 2026 line in the individual spectra were replaced by random 
noise with S/N similar to the uncontaminated regions. We show the 
average spectrum around Zn II A 2026 A in Figure 1 1 . A feature 
was seen in the region with Wo=8±2 mA and log Nzn ii=11.63. 
The feature at -^-1-70 km s^^ may be due to absorption from Mg 
I A 2026.477, and was not included when determining Wq. The 
systems in this sample have an average H I column density of log 
<Nhi>=19.80, giving this sample an average Zn abundance of 
<[Zn/H]>= —0.80 based on this composite spectrum. To increase 
the S/N, we also binned the averaged spectrum to a dispersion of 
0.5 A/pixel ( R'--^4000 ). After binning, a feature was detected at a 
>2(T significance level with Wo=llmA. If the feature is assumed 
to be a true detection, the Zn II column density is log Nzn ii=l 1 .80, 
with a corresponding mean abundance of < [Zn/H] >=— 0.63 for 
these absorbers at 0.7 < Za(,3< 1.5 based on the binned spec- 
trum. In comparison, DLA absorbers have a mean Zn abundance 

of < [Zn/H]>^ 1.0 in this redshift range jKulkami et^l2005l 

l2007h . 

The Nh i -weighted mean metallicity given by: 



Averaged Spectrum WQ=8±2mA 6 v = 30.75 km s 



[(X/H)] = log((X/H)) - log(X/H)0 



where 



((X/H)> = 



ELiN(X)i 
ELiN(H)i 



(2) 



(3) 



has often been used as a quantitative way of estimating the 
amount of metal enr i chmen t in the Universe in a given epoch 
jKulkarni et id] |2005|. l2007h For the absorption systems in this 
sample, the Nh i -weighted mean metallicity based on Fe is 
[<Fe/H>]=— 0.76. Using Zn which is less affected by dust deple- 
tion, we determine the Nh i -weighted mean metallicity for these 




RodiQl Velocity km s 

Figure 11. The unbinned average spectrum of the 10 absorbers in this sam- 
ple near Zn II A 2026 A. The horizontal red line indicates the continuum 
level, and the horizontal dashed green lines show the best fit continuum 
level increased or decreased by ±40% of the RMS variance in the noise. 
The vertical dashed lines show the integration limits. The feature seen was 
detected at > 4(t with Wo=8±2 mA. 



systems at 0.7 < Zai,s< 1.5. to be [<Zn/H>]=— 0.40 assuming the 
upper limits to be detections, and [<Zn/H>]=— 0.70 assuming the 
limits to be zero. 

In this paper we have presented medium-resolution spectra 
of 10 sub-DLAs at Although to date the DLA systems 

have been the preferred tracer of metallicity at high redshift, ev- 
ery absorber with solar or higher metallicity that has been de- 
tected has been a sub-DLA ( Pettini et al. 2000; Khare et al. 200i; 
|PerouxetalJl2006al : IProchaska et alj 12006 : IMeiring et al.l 120070 . 
With this latest sample of sub-DLAs, we have found two more 
systems with high metalhcity, [Zn/H]= -0.05 and [Zn/H]>-l-0.86 
Two other systems in this sample that did not have a Zn II detection, 
did nonetheless have [Fe/H]> —0.30. 

These observations suggest that the largely ignored sub-DLAs 
may provide important clues into the chemical evolution of the 
Universe. Clearly, not all sub-DLAs have solar or higher metal- 
licity, nor would one expect them to if they trace typical galaxies 
which have a range of metallicities. Sub-DLAs do appear, however, 
to have a higher mean metal licity and show faster evolution than do 
DLAs jKulkami et al.l2007h . We continue to find systems with low 
Nh I and high metallicities, but higher S/N is needed to determine 
if the systems with upper limits of [Zn/H]~-0.3 are indeed at that 
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level, or are much less, comparable to DLA abundances. With our 
sample sizes of sub-DLAs now increasing, we can begin to address 
some of the questions relating to the chemical evolution of met- 
als. Future UV spectra would also allow measurements of key lines 
of S III and Fe III which could yield important information about 
ionisation. In forthcoming work, and with an extended sample, we 
will examine differences in the properties of sub-DLAs and DLAs, 
including abundances and kinematical structure, as well as the red- 
shift evolution of these properties in sub-DLAs. 
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APPENDIX - Nh I DETERMINATIONS 

The systems studied in this work all have known Nh i from HST 
spectra. For completeness and convenience to the reader, we pro- 
vide plots of the Voigt profiles of the Ly man-a transition using the 
best fi t values of the column density from lRao, Tumshek. &: NestoJ 
Due to the low resolution and S/N of these UV spectra, only 
one component was used in the fits. We show in figures AI through 
AI O the Voigt profiles corresp onding to the column densities given 
by iRao, Tumshek. fc Nestoa ( 12006) and convolved with a Gaus- 
sian instrumental spread function based on a two pixel resolution 
element, superimposed on the archival data from HST cycle 6 pro- 
gram 6577, cycle 9 program 8569, and cycle 1 1 program 9382. We 
note that the normalization, i.e., the continuum fit that we defin e 
may differ from that adopted bv lRao. Turnshek."fc Nestod ( l200d) . 
For our continuum fits, a polynomial typically of order 5 or less 
or a cubic spline was used, and the absorption line itself was ex- 
cluded from the fitting region. Also over-plotted are profiles with 
H I column densities smaller and larger by 0.15 dex than the best 
fit values. For the Zai,s=0.7377 system in Q1436-0051 in which the 
Nh I value is given as an upper limit, we plot this value (log Nh i 
= 18.8), along with profiles with H I column densities smaller and 
larger by 0.15 dex. 

A horizontal bar in the middle of the profile or directly beneath 
it denotes the region in velocity space where Mg II absorption was 
seen. This spread in Mg II absorption is likely a lower limit ot the 
spread in H I, as some components could have detectable H I with- 
out detectable Mg II. It should be noted that it is possible that Nh i 
may be overestimated when using only one component to fit the 
profiles. Multiple components spread over the width of absorption 
seen in Mg II may decrease the total Nh i needed to fit the profile, 
and thus increase the abundances. Higher resolution UV spectra 
covering higher Lyman series transitions would help to resolve the 
effects of the wide spread of the Mg II components, and discern 
abundance differences between components which is not possible 
with the current lower resolution UV spectra. 
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Table 14. Relative abundances and abundance ratios for the observed sub-DLAs. The first line below the heading gives the solar reference abundance for that column. 



QSO 



labs Log Nh I 



[Zn/H] 



[Fe/H] 



[Zn/Fe] 



[Si/Fe] 



[Ca/Fe] 



[Cr/Fe] 



[Mn/Fe] AlIII/Alir Mg II/Mg T Mgll/Allir Fell/Alli: 



[X/Y] 







-7.37 



-4.53 



-2.85 



+0.07 



-1.13 



-1.82 



-1.97 



Q 1037+0028 

Q1054-0020A 

Q1054-0020B 

Q1215-0034 

Q 1220-0040 

Q1228+1018 

Q 1330-2056 



1.4244 
0.8301 
0.9514 
1.5543 
0.9746 
0.9376 
0.8526 



20.04±0.12 
18.95±0.18 
19.28±0.02 
19.56±0.02 
20.20±0.07 
19.41±0.02 
19.40±0.02 



< -0.63 

< -h0.18 

< -0.21 

< -0.56 

< -1.14 

< -0.37 

< -0.07 



-0.67±0.12 
-0.09±0.18 
-1.09±0.02 
-0.64±0.02 

-1.33±0.07 
-0.30±0.02 
-1.07±0.02 



< +0.04 

< +0.28 

< +0.88 

< +0.08 

< +0.19 

< -0.07 

< +1.00 



0.18+0.04 



< +0.52 



< -0.76 

< -0.04 

< +0.36 

< -0.20 

< -0.14 

< -0.55 

< +0.20 



-0.34+0.05 
-0.09+0.07 

< +0.23 

< -0.79 

< -0.42 

< -0.47 

< -0.22 



< -0.96 



Q1436-0051A 0.7377 20.08+0.11 -0.05+0.12 -0.61+0.11 +0.58+0.04 



-0.98+0.03 < -0.41 +0.04+0.03 



Q1436-0051B 0.9281 <18.8 



>+0.86 



> -0.07 +0.94+0.06 +0.75+0.05 -0.79+0.04 < +0.20 +0.22+0.05 



Q1455-0045 1.0929 20.08+0.06 < -0.80 -0.98+0.06 < +0.18 +0.00+0.10 



-0.35+0.13 -0.51+0.15 < -0.18 



>2.28 
>1.68 

>2.69 
>2.29 
>1.88 

>2.63 
>2.39 



>2.29 
>1.24 
>1.61 

>2.19 

>1.62 

>1.72 
>1.27 



5 
I 

Oh" 

>: 
do 
•S 

I 



"Ratios of column densities 



00 



Sub-DLAs and strong Lyman-limit systems at 1-5. 



19 



Q1037+0028 log N = 20.04±0.15 z„j,= 1.4244 



^ 2x10"'^ - 




Q1054-0020 log N = 18.95±045 z„j,= 0.8301 



Figure Al. 

Q1 054-0020 log N = 19.2a±0.15 z„j,= 0.9514 




Wavelength A 

Figure A3. 



G'17?0-004r. 
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1210 1220 
Wavelength A 

Figure AS. 
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Figure A2. 

Q1215-0034 log N = 19.5b±0.15 z„j,= 1.554.3 



^ 1.0x10^^ - 




Figure A4. 

C'12/8+1018 log N = 19.41 ±0.15 z,i,= 0.93 




1210 '220 
Wavelength A 

Figure A6. 



Figure Al - A6: UV spectra of the systems (clockwise from top left) Za6s=1.4244 in Q1037+0038, z„f,s=0.8301 in Q1054-0020, z„i,s=1.5543 
in Q1215-0034 Zaf,s=0.9376 in Q1228+1018, z<j63=0.9746 in Q1220-0040 and Zn63=0.9514 in Q1054-0020. The dashed green line indicates 
the best fit level of the continuum, and the dashed blue lines are increased or decreased by 10% of the best fit value. Superimposed on the data 
are t heoretical Voigt profiles that hav e been convolved with a gaussian instrumental spread function. The middle profile is the best fit value 
from lRao, Tumshek. fc NestoJ ilOOd). and the upper and lower column densities have been modified by ±0.15 dex from the best fit value. The 
horizontal bar either in the center of the profile or directly beneath the line center denotes the region in velocity space where Mg II absorption 
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Q1330-2056 log N = 19,40±0.15 z„j,= 0,8526 Q1436-005' log N = 20,08±0.15 z,^= 0,737- 

1,5X10""| • • \ • • ' ' \ • • • ' \ ' • • ' \ ' • • • I 3xl0""| 'I'M I I I I 




Wavelength A Wavelength A 

Figure A9. Figure AlO. 

Figure A7 - AlO: Same as Figures 1 1-16, but for the systems (clockwise from top left) 2^,53=0.8526 in Q1330-2056, 2^,53=0.7377 in Q1436- 
0051, Zn5s=1.0929 in Q1455-0045. and z„5s=0.9281 in Q1436-0051. 
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